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The high energy neutrino emission from gamma-ray bursts (GRBs) has been expected in various 
scenarios. In this paper, we study the neutrino emission from early afterglows of GRBs, especially 
under the reverse-forward shock model and late prompt emission model. In the former model, the 
early afterglow emission occurs due to dissipation made by an external shock with the circumburst 
medium (CBM). In the latter model, internal dissipation such as internal shocks produces the 
shallow decay emission in early afterglows. We also discuss implications of recent Swift observations 
for neutrino signals in detail. Future neutrino detectors such as IceCube may detect neutrino signals 
from early afterglows, especially under the late prompt emission model, while the detection would 
be difficult under the reverse-forward shock model. Contribution to the neutrino background from 
the early afterglow emission may be at most comparable to that from the prompt emission unless 
the outflow making the early afterglow emission loads more nonthermal protons, and it may be 
important in the very high energies. Neutrino-detections are inviting because they could provide 
us with not only information on baryon acceleration but also one of the clues to the model of early 
afterglows. Finally, we compare various predictions for the neutrino background from GRBs, which 
are testable by future neutrino-observations. 
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I. INTRODUCTION 

Gamma-ray burst (GRB) is one of the the most ener- 
getic phenomena in the universe and one of the candi- 
dates where both electrons and protons are accelerated 
up to very high energies. If protons are accelerated up 
to very high energies, we can expect high energy neutri- 
nos and gamma-rays that are produced by the photome- 
son production process. We investigated the high energy 
neutrino emission from GRBs under the internal shock 
model in the previous work [l[. Since the prediction by 
Waxman & Bahcall [|J , this kind of neutrino emission has 
been discussed by several authors (see, e.g., [!, 0, [H ) . 

The standard (internal-external) shock model of GRBs 
succeeded in explaining many observations in the pre- 
Swift era (see reviews, e.g., [y, @, H|). Synchrotron ra- 
diation from a reverse shock (RS) and a forward shock 
(FS) usually peaks in the infrared-to-optical or even lower 
bands and ultraviolet-to-x-ray bands, respectively. Some 
of the infrared/optical flashes that can be interpreted 
as the RS emission were observed from GRB 990123 
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IITJIOlS, GRB 021004 [laijand GRB 041219a 
2J, [25|, [20| . These infrared/optical photons can inter- 
act with protons accelerated at the RS and generate high 
energy neutrinos. Under the RS model, Waxman & Bah- 
call predicted high energy neutrino afterglows for the ho- 
mogeneous circumburst medium (CBM) (or interstellar 



medium (ISM)) [27j|. For the wind-like CBM, Dai & Lu 
predicted (10 15 -10 17 ) eV neutrinos Dermer [H,[30| 
and Li et al. [3ll ] considered the neutrino emission from 
the FS, assuming that protons can be accelerated via the 
second-order Fermi acceleration mechanism. Now, large 
neutrino detectors such as IceCube [32[ and KM3Net [33| 
are being constructed. By exploiting these detectors, 
we can test various predictions for high energy neutri- 
nos from GRBs in the near future. Furthermore, ANITA 
[34j and Auger [3|| might also be useful for detecting very 
high energy neutrinos. 

Such possibilities for the neutrino emission have been 
studied under the standard model. However, the recent 
observations by the Swift satellite have shown many un- 
expected behaviors in the early afterglow phase and the 
simplest standard model is confronted with difficulties 
(see e.g., [3(||32l|3l,|3^]). For example, one of them is the 
shallow decaying behavior, which appears after the steep 
decay phase. The flux-decay becomes shallow, when the 
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[40j, |4l| . Although various modifications of the standard 
model have been put forward to explain the observations, 
none of which seems conclusive. 

In addition, the central engine may last much longer 
than the duration of the bursts. The variability of some 
GRB afterglows implies such a prolonged activity of the 
central engine [421 ] . Flares that are found in a significant 
fraction of Swift GRBs require a long duration and/or 
re-activation of the central engine, and one of the lead- 
ing models for explaining flares is the late internal dis- 
sipation model OH]. Furthermore, such late internal 
dissipation may explain the shallow decaying behavior 
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[451] . After the earlier prompt emission, there may be a 
tail of activity of the central engine, producing subshells 
of progressively lower power and bulk Lorentz factor for 
a long time. 

The rapid follow-up by the Swift satellite also allows us 
to observe many early afterglows of GRBs in the optical 
bands by ground-based robotic telescopes such ROTSE. 
However, many ground-based and Swift UVOT observa- 
tions have shown that the majority of bursts have very 
dim or undetectable optical afterglows [46| , Although 
foreground extinction, circumburst absorption and high 
redshift provide explanations for many GRBs, there is 
tentative evidence that the strong RS emission is sup- 
pressed. Theoretically, there are some possible reasons 
for this paucity of previously predicted optically bright 
flashes, if they are ascribed to RS emission. One is the 
absence or weakness of the RS. It is because earlier sim- 
plified estimates can over-predict the optical emission 
when the RS is semi-relativistic [USsI and/or the opti- 
cal emission can be much weaker when the ejecta is highly 
magnetized [4!| . Alternatively, Beloborodov pointed out 
the RS emission can be suppressed enough if the prompt 
emission overlaps with the RS emission [50( ■ On the other 
hand, some authors suggested that the RS emission is 
seen as the shallow decay emission [Ell, [52J. They ar- 
gued that significant energy of such the RS emission may 
be radiated in x-ray bands and the shallow decay phase 
can be explained by the RS running into the ejecta of 
relatively small Lorentz factors. 

Although many models have been proposed, there are 
no conclusive ones for explaining early afterglows. More 
and more observations will be needed in order to dis- 
criminate them. High energy neutrino and gamma-ray 
observations can provide us one of the useful clues to 
such models. In this paper, we focus mainly on the high 
energy neutrino early afterglow emission. We calculate 
neutrino spectra more quantitatively than previous works 
under the various early afterglow models. The method 
of calculation is the same as that used in our previous 
works [H, H3, where we use Geant4 (55[ with experi- 
mental data |56l. \El\ . Future observations of high energy 
neutrinos could provide us with not only information 
on baryon acceleration but also physical conditions of 
GRBs. Combined with observations of electromagnetic 
emissions, they would be useful as one of the clues to 
the early afterglow models, although neutrino telescopes 
that are larger than IceCube will be desirable. 



This paper is structured as follows. In Sec. Ill Al we 
briefly review the reverse-forward shock model. Then, we 
consider contributions from the early and late prompt 
emission in Sec. Ill Bl and III CI In Sec. UH we review 
the high energy neutrino emission process in GRBs. The 
numerical results are shown in Sec. IIVI We describe 
summary and discussions in Sec. [V] 



II. THE MODEL 
A. The Reverse-Forward Shock Model 

The original reverse-forward shock model, which is suc- 
cessful in interpreting many late time afterglow obser- 
vations basically, has been developed (see reviews, e.g., 
[1,0, H, E3, and often accepted on interpreting some 
early infrared/optical flashes such as GRB 990123. Here, 
we briefly review this model which is used for our calcu- 
lations [It], [ll| . The analysis in this subsection is entirely 
within the context of the external shock model. We use 
this model only for the afterglow emission, and assume 
that the prompt emission occurs due to internal dissi- 
pation such as internal shocks. The possibility that in- 
ternal dissipation makes the early afterglow emission is 
described in Sec. Ill CI 

The expanding fireball strikes the surrounding medium 
and will form two shocks: a reverse shock and a forward 
shock. The shocked ambient and ejecta materials are in 
pressure balance and are separated by a contact disconti- 
nuity. In the original standard model, the RS is thought 
to be short-lived, which exists during the initial decel- 
eration of the fireball. After the RS crosses the ejecta, 
the FS continues and the ejecta will transit to self-similar 
expansion which is described by the Blandford-Mckee so- 
lution. 

Each shock compresses the fluid. Let the Lorentz fac- 
tor of the shocked CBM and that of the shocked ejecta 
(measured in the frame of the unshocked CBM) be la- 
beled by r. r is the Lorentz factor of the unshocked 
ejecta (measured in the frame of the shocked ejecta), 
which is written as, 



r 



r 



(i) 



where Tq is the Lorentz factor of the unshocked ejecta 
(measured in the frame of the unshocked CBM) and we 
have assumed ro,r 3> 1. By combining this equation 
with the pressure equality at the contact discontinuity, 
we have, 



r„ 



1 + 2r (n/n cj ) 1/2 



1/2- 



(2) 



Here the number density of the ejecta, rt C j is given by 
n C j = _E e j/47rmpC 2 ro(r A)r 2 , where E e - 3 is the isotropic 
energy of the ejecta and A is the geometrical thickness of 
the ejecta measured in the stellar frame. The circumburst 
number density is written asn = Ar~ k , where k = ex- 
presses a homogeneous CBM (ISM) and k = 2 does a 
wind-like CBM. In the latter case, A = 3 x 10 35 A* cm -1 , 
where A* is the mass-loss rate to wind speed ratio, nor- 
malized to 10 _5 M£ ! yr _1 . A* ~ 1 is a typical value for 
Walf-Rayet stars [58f . 

Now, we define the crossing radius r x , at which the RS 
finishes crossing the ejecta and injection of fresh electrons 
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by the RS ceases. The thickness of the ejecta at the 
crossing time t x is written as, A(r x ) = J^ x dr(f3 — /3rs)> 

where (3 — /3rs — ^-^In 1 ^ 2 /Tln^ 1 / 2 is the difference 
between the speed of the unshocked ejecta in the stellar 
frame and that of the RS. This expression is valid for 
a wide range of the ratio Tgn/n C j, from Newtonian to 
relativistic RS cases [l?], [H|. If the ejecta thickness is 
smaller than the spreading radius which is given by r s ss 
2TqcT, the ejecta thickness at the crossing time A x is 
expressed as A x « Ao = cT. This case is called a thick 
ejecta case in the usual terminology [59(. Here, T is the 
duration of the GRB ejecta in the stellar frame. If the 
ejecta thickness is larger than r s /2rQ, we have A x 
r/2T 2 . This case is called a thin ejecta case. We can 
write the crossing radius in the thick ejecta regime as 
follows, 



0.24 x 10 17 (^p^) 1/4 cm (fc=0) 
0.52 x 10 15 (%^) 1/2 cm (k=2) 



(3) 



In the thin ejecta regime, we have the following expres- 
sions, 



The magnetic fields in each post-shock fluid are param- 
eterized by fractions e B and e r B of the post-shock internal 
energy density stored in it, respectively. They are written 



/ 



Bl = 



^32TTe B nm p c 2 (T x - 1)(T X +3/4) (5a) 
^2ire B n e - 3 m p c 2 (T' x - l)(r' x +3/4) (5b) 



The injection Lorentz factor of electrons at the crossing 
time is, 
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where g(p) is given by g(p) — (p — 2)/(jp — 1) and f e is a 
number fraction of the shocked electrons that are injected 
into the acceleration process. The cooling Lorentz factor 
of electrons is obtained by equating the cooling time and 
dynamical time td yn = r x /T x c, which is written as, 



67TTOpC 2 r> 



7e,c = 



( 7T J B 2 x r x (l + r)' 



(7) 



Bo 



0.44 x 10 17 ( 
L0 * 1014 (xfcj 



1/3 



cm (fc=0) 
cm (fc=2) 



(4) 



The essential length scales for the thick and thin ejecta 
cases can be expressed as r& = ^^Ao 1 ^ 4 and rr = 
ffo~ 2 / 3 , respectively. Here, I is the Sedov length. For 

1/3 

k = 0, we have I = (-E G j/(47r/3)nm p c 2 ) . Using 
these length scales, we can write the crossing radius 
r x ~ max(rA,rr). It is often convenient to introduce 
£ = {r T /r A f = {l/A ) 1/2 T- 4 / 3 . Note that, if £ > 1, the 
ejecta is essentially in the thin ejecta regime. If ^ < 1, it 
is essentially in the thick ejecta regime. 

The crossing time is calculated by ( x = 
Jq X dr(l/2T 2 c). For a thick ejecta case with k — 0, we 
(e \ 1/3 

can obtain t x fa 17 s (^ „ r T 53 J ■ For a thick ejecta 

case with k = 2, i x « 2.9 s ( a^t^ ) • ^ or e j ec * a 
cases with k = and k = 2, we have t x « 0.72 s x T. 
Up to an order-unity factor, the crossing time can be 
approximated as t x ~ max[T, tr}- Here, the crossing 
time is regarded as the time when afterglows begin, i.e., 
the self-similar behavior starts. The RS emission rises 
after the ejecta starts to sweep up matter, but it does 
not peak until time ~ t x , marking the beginning of the 
self-similar regime. Most of the RS emission (hence the 
neutrino emission from the RS) is released at ~ t x . 
Therefore, we can expect that it is sufficient to consider 
the behavior at t x for our purpose of evaluating the 
neutrino energy fluence, although more sophisticated 
calculations including the time dependent evolution are 
possible [6lj . 



where Y is the Compton parameter, which can be calcu- 
lated from the electron distribution (see e.g., [60|). 

The injection energy e m = hv m and cooling en- 
ergy e c = hv c in the comoving frame are calculated 
from the corresponding electron Lorentz factors by 



hv 



2eB y 



In this paper, we con- 



usmg £ , 

sider the four cases; the thick ejecta colliding into 
the ISM, the thin ejecta colliding into the ISM, the 
thick ejecta colliding into the wind-like CBM, and 
the thin ejecta colliding into the wind-like CBM. For 
example, let us consider the case of the thin ejecta 
colliding into the ISM. We can obtain two character- 
istic frequencies in the observer frame (v oh ~ r x ^) 
as follows. For the FS emission, we have z/jj 3 ^ ~ 

2.9 x lO^Hzs 2 ^^,) 172 ^^)^//)- 2 ^^ 1 ^ and 

^ » 4.5xl0 16 Hz(4 i _ 2 )- 3/2 n - 1 S cj , 1 5 / 3 2 T- 1 / 2 (l + Y)-\ 
For the RS emission, we have v^r ~ 1-0 x 
lO^Hzg 2 ^) 1 / 2 ^) 2 ^)" 2 ^ - l) 2 n 3/8 ^ 5 8 3 T- 1 /8, 

< b r « 4.5 x W™R Z (e B )- 3/2 n^E^% 2 T-^(l + Y)- 2 . 
Basically, the observed cooling frequencies of RS and 
FS emissions are equal if microphysical parameters are 
similar, while the observed injection frequency of the RS 
emission is smaller than that of the FS emission by a 
factor of T 2 . 

We also consider the synchrotron self-absorption pro- 
cess. The optical thickness for synchrotron self- 
absorption can be approximated by (see, e.g., [60|). 



T sa (e) 



be 



A-KT 2 x B xl l n 



fN xJ (£/e " r " {£<£r ' 

J e l\ e X < (p + 4) 
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(e™ < e) 
(8) 



4 



where 



n: = 



47T , 

Ar 

3-fc 

T m p c 2 



3-fc 



(9a) 
(9b) 



are the number of electrons energized by the FS and RS, 
respectively. 7 e ,„ = min(7e iC , j e ,m), £ n = min(e m ,£ c ), 
and we have assumed 7 ejn > 1. The self- absorption en- 
ergy is determined by T sa (e sa ) — 1. Hence, we can obtain, 



T,«*(e = e n ) (£<£") 
r sa ^(e = s n ) (e n <e) 



(10) 



Generally, the observed self-absorption frequency in the 
RS is typically higher than that in the FS. 

The comoving specific luminosity per unit energy at 
the injection energy or cooling energy is approximated 
by 



= 1 f c N c ^ e3B : 
2ith rripC 2 



'pi 
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where 



is an order-unity coefficient calculated by 



Wijers & Galama [62|. In this paper we use 
<j) p ~ 0.6. For example, let us consider the 
case of the thick ejecta colliding into the ISM. 
The observed peak luminosity per unit energy from 
the FS emission is written as Ll _„ = 1.3 x 

10 58 s- 1 //^ cj ,53(4_ 2 ) 1/ ^o /2 T 1 / 2 ( ( /) p /0.6), and that 
from the RS emission is expressed as L r e max = 4.9 x 

10 62 s" 1 /e r B£(ea 1 / 2 n;/ 4 r o -Lr : 1 / 4 (^/0.6). The 
observed peak flux from the RS is typically larger by 
a factor of T x than that from the FS. 

From the above equations, we can obtain photon spec- 
tra. We have for e sa < e m < e c (the slow cooling regime), 



dn 
de 



{e sa /e m y^{e/e sa ) 1 



( e min ^ e ^ e sa ) 
{e sa < e 



<e^e c ) 

(e c /s m y E ^~(e/e c y E ¥ ( e c < £ ^ e max ), 

(12) 

and for e sa < e c < e rn (the fast cooling regime), 



dn 

ds 



{e sa /e c y T ye/e sa y 
(e/eT 1 3 

(e m /e c y^(e/e m y E ^ 



Similarly we can obtain for e m < e & 



(e min ^ e S £ sa ) 
(e sa <e^e c ) 
(e c < e ^ e m ) 

(e m < e ^ e max ). 

(13) 

<£ c , 



dn 
de 



{e sa /e m y E ^ 1 {e m /e sa )^{e/e 7n ) 1 

{e sa /e m y E ^ 1 {e/e sa )i 
p+i 

{e/e m y — 

P+l p+2 

y (s c /e m y—(s/e c y 2 , 



(14) 



and for e c < e sa < e m (where we do not consider the 
inhomogeneity [63j|). 



dn 
de 



where 



{e sa /e c y^(e c /e sa )^{e/e c ) 1 
{e sa /e c y^{e/e sa °y 

(£,£Cyi 3 

(e m /e c y^(e/e m y E ^ , 



(15) 



47rririce r; 



(16) 



In order to demonstrate neutrino spectra from the RS, 
we adopt following parameter sets. 

ISM-tc: r = 10 2 - 5 , £ C j = 4xl0 53 ergs, A = 4.5 xlO 11 
cm, n = 5cm- 3 , e r B = 0.01, e r e = 1/4, f r e = 1 andp = 2.4. 
This is one of the parameter sets representing the thick 
ejecta that collides into the ISM. 

ISM-tn: r = 10 2 , E ci = 4 x 10 52 ergs, A = 4.5 x 10 11 



cm, n = 0.5 cm 



0.01, e r e = 1/4, fl = 1 and 



p = 2.4. This is one of the parameter sets representing 
the thin ejecta that collides into the ISM. 

ISM-e: T = 300, E c] = 10 55 ergs, A = 1.05 x 10 12 
cm, n = lcm- 3 , e r B = 0.001, e r e = 0.04, fl = 1 and 
p = 2.5. This is one of the parameter sets representing 
the thin ejecta that collides into the ISM [l7l |. 

ISM-eb: T = 300, E c] = 10 54 ergs, A = 1.05 x 10 12 
cm, n = 1 cm~ 3 , e r B = 0.2, e r e — 0.1, /J = 1 and p = 2.5. 
This is one of the parameter sets representing the thick 
ejecta that collides into the ISM [l6|. Note that the mod- 
els ISM-e and ISM-eb express cases where the ejecta is 
energetic E Ci > 10 54 ergs. GRB 990123 is thought as 
such an energetic event. (Note that, if the ejecta car- 
ries the magnetic field directly from the central engine, 
e r B can be significantly larger than e B . A strong reverse 
shock is expected when the magnetic field is radiation- 
ally important but not yet dynamically important. GRB 
990123 is thought to be one of such events [i"6|.) 

WIND-tc: T = lO 2 ' 5 , E ci = 4 x 10 53 ergs, A = 
4.5 x 10 11 cm, A* = 1, e r B = 0.01, e r e = 1/4, f r e = 1 and 
p = 2.4. This is one of the parameter sets representing 
the thick ejecta that collides into the wind-like CBM. 

WIND-tn: r = 10 2 , £ cj = 4 x 10 52 ergs, A = 4.5 x 
10 11 cm, A* = 0.01, e r B = 0.01, e r e = 1/4 and p = 2.4. 
This is one of the parameter sets representing the thin 
ejecta that collides into the wind-like CBM. 

So far, we have assumed f e — 1. However, this might 
not be true (64|. We can consider cases with / e < 1. 
If a small fraction of the electron population can have a 
significant part of the dissipated e nerg y in the RS, the 
RS emission can appear in x-rays [5l| . The RS emis- 
sion, which has usually been used for explanation of in- 
frared/optical flashes that are observed for only a small 
fraction of GRBs, might play an important role in the 
early afterglow phase. Recently, some authors proposed 
that the mysterious shallow decay emission can be ex- 
plained by the RS emission [5l|, [52| • In their models, the 
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plateau shape can be achieved by requiring the appro- 
priate distribution of Lorentz factors of the ejecta. It is 
assumed that Lorentz factors of the material, which is 
ejected during the last stages of source activity, decrease 
to small values of Lorentz factors. The head of the ejecta 
has larger Lorentz factors, while the tail of the ejecta has 
smaller Lorentz factors. In such long-lived RS is 

possible when it propagates into the stratified ejecta with 
decreasing Lorentz factors. In addition, this model also 
requires that the forward shock emission can be negligible 
at least in the early afterglow phase. The suppression of 
the FS emission might occur because the magnetic fields 
are too weak in the external medium and not sufficiently 
amplified [HI], and/or the first-order Fermi acceleration 
mechanism is not so efficient. To reproduce the shallow 
decaying behavior, we need the detailed numerical mod- 
eling. This is beyond scope of this paper because our 
goal is not to explain the shallow decay emission. For 
our purpose to estimate neutrino fluxes, it will be suf- 
ficient to consider the RS emission by the head of the 
ejecta. Assuming that the head of the ejecta carries the 
energy E^ = (l/3)E e j, we adopt the following parameter 
sets. 

ISM-s: T = 10 2 ' 5 , = | x 10 53 ergs, A = 4.5 x 10 11 
cm, n = 5 cm- 3 , e r B = 1/4, e r e = 1/4, f r e = 0.025 and 
p = 2.4. This is one of the parameter sets representing 
the thick ejecta that collides into the ISM. 

WIND-s: r = 10 2 - 5 , ££j = | X 10 53 ergs, A = 4.5 x 
10 11 cm, A* = 0.1, e r B = 1/4, e r e = 1/4, f r e = 0.025 and 
p = 2.4. This is one of the parameter sets representing 
the thick ejecta that collides into the wind-like CBM. 

In order to predict neutrino fluxes, we have to set the 
amount of accelerated protons which no one knows from 
the first principle. In this paper, we just assume that the 
moderately efficient acceleration occurs and take e acc = 
1/4, where e acc = ( p (l — €b — £e) and £ p is acceleration 
efficiency. 

B. The Overlapping of Prompt Emission with the 
Shocked Region 

In the previous subsection, we have considered the RS 
emission under the original RS model, which typically 
predicts infrared/optical flashes. As previously noted, 
there is tentative evidence of the lack of infrared/optical 
flashes jigf . Several possible reasons have been sug gested. 
First, the ejecta may be strongly magnetized [49]. Then, 
the hydrodynamical shock can become weak or there is 
no RS. Second, especially in the thin ejecta, the RS emis- 
sion can be more suppressed than the earlier simplest es- 
timations [13, 13 . Such a semi-relativistic RS may give 
a peak flux below one give by a FS. In addition, a pair- 
rich RS may be common (48[. Third, in the thick ejecta, 
The RS emission can be suppressed because RS electrons 
are rapidly cooled due to Compton scattering by photons 
from the prompt emission [5fl |. 

Now, we consider the third possibility, i.e., overlapping 



of the prompt emission with the shocked region. This 
overlapping will occur for the thick ejecta case, which 
can be expected for long GRBs we consider throughout 
this paper. The prompt emission which occurs due to 
internal dissipation can provide additional target photons 
for accelerated protons, so that more neutrinos can be 
produced via the photomeson production process. This 
possibility was suggested by Fan et al. [I(| for the thick 
ejecta colliding into the wind-like CBM. Here, we also 
study such possibilities including cases of the ejecta that 
collides with ISM in more detail. We assume E™° = 
10 53 ergs as the isotropic prompt emission energy. The 
averaged photon energy density within the ejecta at the 
crossing radius is given by, 

7 4^r 2 A r 2 v ' 

A spectrum of the prompt emission is well approximated 
by a broken power-law spectrum, which is, 

dn f e' a (for e min < e < e b ) . . 

de^X e-P (for e b < e < e max ). 1 j 

The observed break energy is ~ 250 keV, which corre- 
sponds to the break energy in the comoving frame, e b ~ 
a few keV. Hence, we set e b — 1 keV and spectral indices 
to a = 1, (3 = 2.2 similarly to our previous work [l|. 
We take the minimum energy as 1 eV and the maximum 
energy as 10 MeV in the comoving frame. 

If such a photon flow due to the prompt emission 
cools down high energy electrons sufficiently, strong in- 
frared/optical flashes that are expected in the RS model 
can be suppressed In the model ISM-tc which rep- 
resents the thick ejecta case, the cooling Lorentz factor 
of electrons 7 EiC is significantly lowered due to Compton 
scattering by prompt photons. It leads to that the RS 
emission is not in the slow cooling regime but in the fast 
cooling regime. 

In this paper, we consider the four parameter sets; 
ISM-tc2, ISM-s2, WIND-tc2 and WIND-s2. Model pa- 
rameters are the same as those in ISM-tc, ISM-s, WIND- 
tc and WIND-s, respectively. But the overlapping effect 
due to the prompt emission is included. A target photon 
spectrum is given by a superposition of a prompt spec- 
trum and a RS spectrum modified by the overlapping ef- 
fect. Note that the high energy gamma-ray emission can 
be expected in this model [50|, |66[. These up-scattered 
photons can also contribute to the photomeson produc- 
tion process, but we can neglect such a population be- 
cause of a smaller number of these high energy photons. 

C. The Late Prompt Emission Model 

The late internal dissipation may last longer than the 
duration of earlier internal dissipation that makes the 
prompt emission. For example, some of the flares are 
likely to be attributed to the late internal dissipations 
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|43l l44j • The flares typically happen hundreds of seconds 
after the trigger or earlier. In some cases, they occur 
around a day after the main burst. The amplitudes of 
the flares are usually larger than the underlying afterglow 
component by a factor of several, but can be much larger. 
Recent analyses show that the averaged radiation energy 
of flares is approximately ~ 1/10 of that of the prompt 
emission [44j . 

One of the leading models for such late internal dissi- 
pation is the late internal shock model (see, e.g., [13] )• 
Lorentz factors of ejected subshells will be highly vari- 
able. If r s ~ 10 and Tf ~ 100 are typical Lorentz 
factors of the slow and fast subshells respectively the 
Lorentz factor of the merged subshell can be expressed 
as Lo ~ VTfTs ~ 30. The Lorentz factor of the internal 
shocks can be estimated by L s h ~ (L//L s + L s /L/)/2 ~ a 
few. The typical collision radius is given by, r ps 2T^c5t 
10 15 - 3 cm (L /15) 2 [<5i/150(l + z) s]. The internal shocks 
are expected to be mildly relativistic shocks. If protons 
are accelerated efficiently in these shocks, neutrinos can 
appear through the photomeson production process. In 
our previous work [53j, we predicted such high energy 
neutrino flashes from flares of GRBs. 

On the other hand, the x-ray emission in the shallow 
decay phase is often attributed to the external shock 
emission. However, Ghisellini et al. (4|| recently sug- 
gested that this plateau phase for the x-ray emission may 
be due to the late prompt emission. In their model, this 
late prompt emission can be due to the same internal 
dissipation process as that for the early prompt emission 
but by subshells created at late times with smaller r 
and much lower power. The radiation can be produced 
at distances relatively close to the central engine (even 
less than r ~ 10 13 ~ 14 cm), in a different region where the 
subshells interact with the CBM |45j . 

Although this model can explain the chromatic behav- 
ior in early afterglows, it has not enabled us to explain 
the closure relations, which are expected in forward shock 
models and satisfied for the normal decay segment fol- 
lowing the shallow decay segment. In spite of such a 
defect, there are indeed a couple of bursts that are likely 
to show the late prompt emission, marked by a sharp de- 
cay following an extended plateau with flickering. Such 
a strikingbehavior is seen in some bursts such as GRB 
070110 [6j|. 

If the shallow decay emission is attributed to the late 
prompt emission and late internal dissipation is due to 
late internal shocks or other models that allow a signif- 
icant fraction of baryons to be accelerated up to suffi- 
ciently high energies, we can also expect high energy neu- 
trino signals coincident with the shallow decay emission 
like the case of flares. The neutrino emission due to late 
internal dissipation itself would decay during the shallow 
and normal decay phases after the late prompt emission 
starts. Note that, physical conditions may be similar to 
those expected in flares and flares may be produced by 
a late shell, moving with a somewhat larger Lorentz fac- 
tor than the shells created just earlier [45|]. Therefore, 



we can use the same framework as that used in Murase 
& Nagataki [53[ in order to evaluate high energy neu- 
trino fluxes associated with the early x-ray emission in 
the shallow decay phase under the late prompt emission 
model. In this paper, let us refer neutrinos from both 
flares and the late prompt emission as neutrinos from 
the late prompt emission. 

The photon energy density is given by, 



4tTV 2 Z 



(19) 



where I is the width of subshells and E™° h is the radi- 
ated energy from each subshell. I is typically given by 
I = r/To, although it can be a smaller value. We assume 
the relatively small Lorentz factor, Lo ~ a fewxlO. The 
emitted energy in the shallow decay phase is typically 
~ 1/10 of that of the prompt emission [69j]. Hence, we 
adopt the total isotropic radiation energy of the x-ray 
emission in the shallow decay phase, E^p = NE^° h = 
10 52 ergs, where N is a number of collisions between 
ejected subshells. The shallow decay emission typically 
occur at T ~ 10 3 ~ 4 s, and a collision radius can be 
r ~ 10 13 ~ 16 cm under the late prompt emission model 
(and note that the expected variability time scale under 
the late internal shock model can be estimated by using 
St fa (1 + z)r/2Tlc). Hence, we take N ~ (a few-100). 

For numerical calculations, we assume a broken power- 
law spectrum similarly to the cases of the prompt emis- 
sion, and we use Eq. (18). But we take e b = (10 — 100) 
eV. We also set e min = 0.1 eV and e max = 1 MeV. 

Magnetic energy density is expressed as Ub = £bU 7 
where £g ~ £s/e e . In this paper, we set £g = 1. 
Nonthermal proton energy density is parameterized as 
Up = £acc£/ 7 , where £ acc ps e acc /e e = C, p {l-e B -e e ) / 'e e is a 
nonthermal baryon loading factor, where £ p is the proton 
acceleration efficiency. The present acceleration theory 
cannot give this value from the first principle. Here, we 
assume that protons can be accelerated efficiently as is 
usually expected in supernova remnants and just adopt 
£acc = 10 as a fiducial value. If the proton distribution 
has drip/dEp oc e~ 2 and the minimum energy of protons 
is a fewxiripC 2 , this fiducial value £ acc = 10 means that 
the energy density of nonthermal protons with energy e p , 

er p(de E ) * s comparable to the radiation energy density 

of photons with energy e b , ( £& ) (df) e _ e &' A similar as- 
sumption is often adopted in previous works P, 0, d, [f| . 
If the nonthermal baryon loading factor could be larger, 
we expect higher neutrino fluxes, but no one obtains this 
value from the first principle. On the other hand, there 
is a hypothesis that UHECRs come from GRBs based 
on the internal shock model that will cause the usual 
prompt emission ,2], which typically requires sufficiently 
large nonthermal baryon loading factors £ acc ~ (50 — 100) 
[l| (although statements depend on the evaluation of the 
local GRB rate which has some uncertainties [?0j)- Mo- 
tivated by this hypothesis, we also use £ acc = 50 as the 
optimistic value in this paper. 
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We adopt the following parameter sets to estimate the 
neutrino flux from GRBs in the late prompt emission 
model. 

LPO: E™° h = 10 51 - 2 ergs, T = 15 r = 10 15 3 cm, £ B = 
1 and e b = 10 eV. 

LP1: E^° h = 10 51 ergs, r = 10 r = 10 15 cm, £ B = 1 
and e b = 100 eV. 

LP2: E^° h = 10 50 ergs, r = 10 r = 10 14 cm, £ B = 1 
and e b = 100 eV. 

Note that the model LPO corresponds to the model 
FUV-ray flare (A) used in Murase & Nagataki 



III. NEUTRINO PRODUCTION IN GRBS 

We have assumed that the early afterglow emission 
comes from high energy electrons that are accelerated via 
some dissipation process such as shock dissipation. Pro- 
tons also may be accelerated in internal and/or external 
shocks. If the first-order Fermi acceleration mechanism 
is realized, a proton spectrum can be written as, 



dn p 



de 



f dnp_ 
P \ b P de p 



(20) 



where p is the spectral index. Its value is typically p«2 
in the non-relativistic shock diffusive acceleration the- 
ory. In the ultra-relativistic shock limit, p « 2.2 is ob- 
tained assuming the isotropic diffusion in the downstream 
fn\ . [72[. In this paper, we adopt p — 2 for the proton 
spectrum. U p is the energy density of nonthermal pro- 
tons, which is given by U p = £ a ccKy in the late prompt 
emission model while by U p = e acc (E e j/4:Trr' x r x A x ). 
(Note that the neutrino background can be calculated, 
given the GRB rate history. For this purpose, we also 
use the different normalization based on the observed 
UHECR flux. See the end of this subsection. ) The mini- 
mum energy of protons would be £™ m ~ a fewxr re im p c 2 , 
although the exact value is unknown. Here, r re i is the 
relative Lorentz factor, which is r rc i = T' for the RS 
or r ro i = r s h for the internal shocks. Here, we take 
£ mm _ 20 GeV, although the accurate value is irrele- 
vant for the resulting spectra. For the FS, the proton 
energy in the observer frame is given by i?™ m ~ a few 
x Tm p c 2 . 

The maximum energy of cosmic-ray nuclei is deter- 
mined by several criteria. One of the necessary conditions 
is obtained by comparing the Larmor radius of cosmic- 
ray nuclei with the size of the acceleration region [73|, [74j . 
In cases we consider, this criterion corresponds to com- 
paring the acceleration time scale t acc with the dynamical 
time scale tdyn- In addition, the maximum energy of par- 
ticles is limited by various cooling processes and diffusive 
losses of them. Other criteria are obtained by compar- 
ing the acceleration time scale with various cooling time 
scales (the synchrotron cooling time scale t syn , adiabatic 
cooling time scale t a d and so on) and with the escape 



time scale due to particle diffusion t esc 0, llH . For 
cosmic-ray nuclei to be accelerated, all the above criteria 
should be satisfied. As an example, let us consider the 
first criterion and estimate the possible maximum energy 
of protons accelerated in the FS for the thick ejecta col- 
liding into ISM. We have [zl, 



E™ ax » ZeT x B 1SM r> 



3 ISM p 3/8 -3/8^1/8 



5.1xlO iJ eVZSi^ 53 n 



T?°, (21) 



where B 1SM is the strength of the upstream magnetic 
field. From the above equation, we can expect that pro- 
tons cannot be accelerated up to ultra high energies at 
the FS by the first-order Fermi acceleration mechanism 
[HI, [lH . Although the second-order Fermi acceleration 
might allow protons to be accelerated up to ultra high 
energies [76| . we consider only the first-order Fermi ac- 
celeration mechanism for the FS in this paper. From Eq. 
(21), we can see that neutrinos, which are produced by 
the photomeson production process in the FS, will be 
negligible. Corresponding to Eq. (21), the maximum 
energy of protons accelerated in the RS is written as, 



E™* « ZeB r x r x 

« 2.0xl0 21 eVZe 1 J 2 E 3 c l 8 53 n 1 /S T 1 1 / 



(22) 



Therefore, protons can be accelerated up to ultra high 
energies at the RS by the first-order acceleration mech- 
anism if other cooling processes are not important. On 
the other hand, in the late prompt emission model, the 
UHECR production is typically impossible. Note that 
the photomeson cooling process can be important in very 
high energies in this model [53|- We treat various cool- 
ing time scales properly in our numerical results that are 
shown later (see Appendix A). 

Sufficiently accelerated protons can interact with tar- 
get photons via the photomeson production process and 
produce high energy pion and muons. Neutrinos are pro- 
duced via the decay of 7r ± — * /i ± +^ /J (P A1 ) — * e ± + iy e (D e ) + 
v^+Vix. The neutrino production efficiency is represented 
by the photomeson production efficiency f pl = idyn/Vy 
(see Appendix A). For example, let us consider the thick 
ejecta colliding into the ISM. In such cases, we can ap- 
proximately obtain 0, H3, HH , 



fjyy ~ 0.088 



J 6.48 



f (E p /E b f 1 (E p < El) 
r x ,i 6 T Xt2 2 El QeY \ (E p /E b p f 1 (E b p < E p ) 

(23) 

where E b is the observed break energy which is either of 
£ ob or e ob or e ob' an d Lb is the observed luminosity at 
the break energy. E b ~ Q .M ^mp^Y 2 , / E h is the proton 
break energy, where, £a is around 0.3 GeV. For example, 
E b = 10 eV and T x = 100 lead to E b p - 10 20 eV. In Eq. 
(23) we do not include the effect of multi-pion production 
which is not important for the RS emission. 

Next, let us consider the overlapping of the prompt 
emission with the shocked region. We can obtain (includ- 
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ing the effect of multi-pion production which is moder- 
ately important in very high energies) [53l ] . 



in 



1.7 x 10' 



TV* TTiiso 

1 -^7,53 



_ (E p /E b p 
r^i e A 0tl2 Tl 25 E b loakcV \ (E p /E b p ) c 



fA/3" 1 



(24) 

From Eq. (24), we expect that photomeson production 
is not so efficient when the CBM is the ISM. If the CBM 
is wind-like, it becomes more efficient (66j . 

Finally, let us consider the late prompt emission model. 
The physical conditions will be similar to those of flares. 
From Eq. (3) in Murase & Nagataki [53[, we have, 



in 



5.2 



TplSO 

-^7,811,50.5 

r 2 E b 

' 14.5-^1 kcV 



(E P /E b f 



(E P /E, 



(25) 



From Eq. (25), we can expect that almost all the pro- 
tons that are accelerated to sufficiently high energies will 
be depleted due to photomeson production. The above 
approximate evaluation by using Eqs. (23-25) is good 
agreement with numerical results. 

We treat photomeson production in detail and calcu- 
late neutrino spectra numerically. Generally, cooling pro- 
cesses of pions and muons are important. Of course, our 
numerical calculations take into account them, but the 
analytical consideration is convenient. Hence, we briefly 
review the neutrino emission process here. More detailed 
discussions can be found in e.g., Rachen & Meszaros [73l |. 

One of the important cooling processes is synchrotron 
cooling of pions and muons. The synchrotron break for 
neutrinos from pions/muons is determined by t^^^n = 
7w / ' pJ-K i fj,i where t^i^ is the mean life times of pions and 
muons. For neutrinos from pions and muons, we obtain, 



El 



E^ 



67rmijc 5 



4 V cnmlB 2 ^ 



67rm^c 5 



3 r V tTTinrU-T,, 



(26a) 



(26b) 



Above the synchrotron break, a neutrino spectrum is 
suppressed by iw^syn/iw/!- The adiabatic break for 
neutrinos from pions/muons is determined by t^/^ ad = 

Ttt/^t/m- We have ' 



El 



-I — m T c 

4 Tir 



Er « ir^m.c 2 . 

3 T u 



(27a) 
(27b) 



Above the adiabatic break, a neutrino spectrum is sup- 
pressed by ^/^ad/ivr/V 

Now, we can obtain approximate neutrino spectra of 
(vp + Dp). For E b v < E s u {< E%), we obtain, 



E: 



,dN v 



1 



,dN n 



{EjElf- 1 



{Et/ElT \E V /Ety 



(28) 



For E h v < El < E s v , we have, 



E 



,dN v 
* dE v 



1, F2 dNp 
4 M p dE p 

(Ev/Etf- 1 
(Ev/Elf- 1 

{El/Elf \E V /El) a - 2 
(El/E b u f \El/El) a -\E v /El 



(29) 



,c*-3 



Note that the above two approximate expressions are 
valid as long as / P7 < 1. Although more general ex- 
pressions can be obtained, we do not show them here. 
We can obtain general neutrino more quantitatively by 
grace of numerical calculations. 

It is important to see a lot of GRB events, because the 
detection of neutrinos signals from one GRB is generally 
not easy We can calculate the neutrino background 

assuming the GRB rate history (see Appendix B). In this 
paper, we adopt two different normalizations of the pro- 
ton flux. One is the normalization by the observed GRB 
rate and isotropic energy with the fixed acceleration effi- 
ciency. As noted before, we adopt e acc = 1/4 for the RS 
model and £ acc = 10 for the late prompt emission model 
as fiducial values, respectively. The other normalization 
is based on the hypothesis that observed UHECRs (whose 
flux is E^dNp/dEpiz = 0) w 0.65 x 10 44 ergs Mpc^^r" 1 ) 
come from GRBs. When we adopt this normalization, we 
simply change the normalization factor just for simplic- 
ity. For more detailed calculations, proton spectra after 
the propagation in the universe should be fitted to the 
observed UHECR flux, and we have to take care of af- 
fections to target photon spectra by the change of the 
released energy per GRB according to the assumed local 
GRB rate. In the RS model we expect that UHECR pro- 
duction is possible, so that neutrino fluxes from the RS 
should generally satisfy the Waxman & Bahcall (WB) 
bound [?3 because of f pl < 1 typically. For the late 
prompt emission model, we adopt £ acc = 50 as an opti- 
mistic value, as is noted in the previous section. Note 
that neutrino fluxes in this model do not have to satisfy 
the WB bound generally (rather, it should be constrained 
by the more general bound [78[). This is because protons 
are usually not accelerated up to ultra high energies and 
we can expect f pl > 1 in this model. 



IV. NUMERICAL RESULTS 

A. Neutrino Spectra from GRB Early Afterglows 

We need to evaluate proton cooling time scales which 
include the photohadronic cooling calculated by Geant4. 
In Fig. 1, we show the acceleration time scale and pro- 
ton cooling time scales for the model ISM-e which rep- 
resents the case of energetic ejecta that is expected in 
GRB 990123. In this case, the most important cooling 
process is adiabatic cooling. The maximum energy is de- 
termined by Eq. (22). For the ISM cases such as the 
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FIG. 1: The acceleration time scale and various cooling time 
scales of protons for the model ISM-e. The energy scale is 
measured in the comoving frame. 



models ISM-tc and ISM-tn, the adiabatic time scale is 
typically the most important. For the wind-like CBM 
cases, synchrotron cooling can be more important. On 
the other hand, photohadronic cooling can be dominant 
cooling process in the late prompt emission model, as 
shown in Murase & Nagataki [53[. Synchrotron cooling 
and adiabatic cooling can also be important. 

In the case of the prompt emission, the effect of pion- 
multiplicity and proton-inelasticity is moderately impor- 
tant in the sufficiently high energies. This comes from 
the fact that more and more target photons can interact 
with protons via double- and multi-pion production, as 
the energy of incident protons becomes high. As a re- 
sult, the photohadronic cooling time scale increases with 
the energy more than expected in the A-resonance ap- 
proximation Similar things can be applied when we 
consider the overlapping effect of the prompt emission 
and the late prompt emission model. This effect is more 
important for a flatter photon spectrum as demonstrated 
in the previous paper. On the other hand, the RS emis- 
sion (a ~ (1.5 - 1.7) (fore < max[e m ,e c ]) and [3 ~ 2.2 
(fore > max[e m , e c ])) has a steeper spectrum than that 
of the prompt emission (a = 1 and /3 = 2,2). There- 
fore, the multiplicity effect is less important in the case 
of the RS emission and the A-resonance approximation 
is a good approximation. In fact, Eq. (23) is good agree- 
ment with numerical results. 

By executing numerical calculations, we can obtain 
neutrino spectra (see Appendix A) . The adopted param- 
eter sets are summarized in Tables I and II. First, we 
show the neutrino energy fluence from one GRB event. 
The neutrino energy fluence can be evaluated by, 



Ei<b v = 



(1 



z) 2 <m 

de v 



(30) 
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TABLE I: The adopted parameter sets in the reverse-forward 
shock model. Note that the models, where the overlapping 
effect due to the prompt emission is included, are referred 
as ISM-tc2, ISM-s2, WIND-tc2 and WIND-s2. The isotropic 
prompt emission energy _B^ so = 10 53 ergs is assumed in these 
four parameter sets. 



Model 


To 


E^ sh [ergs] 


r [cm] 




E b ' [keV] 


LP0 


15 


10 bi.x 
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1 


0.15 


LP1 


to 


10 5i 


10 lb 


1 


1 


LP2 


10 


10 bu 


to 14 


1 


1 



TABLE II: The adopted parameter sets in the late prompt 
emission model. 



given cosmology. Throughout the paper, we assume the 
ACDM universe with ft m = 0.3, Q k = 0, = 0.7 and 
H = 7Ikms- 1 Mpc~ 1 . 

In Fig. 2, we show neutrino energy fluences for the 
model ISM-e. Even for prompt neutrino bursts predicted 
under the internal shock model, it is difficult to see high 
energy signals from one GRB event by IceCube, which 
typically requires E^c/jl, > I0~ 4 ergcm~ 2 for 0.1 PeV neu- 
trinos |4j]. Therefore, it is much more difficult to see 
neutrino afterglows under the original RS model because 
the smaller number of neutrinos with higher energies is 
expected in this model. 

We comment on the effect of double- and multi- pion 
production to neutrino spectra. In the case of the prompt 
emission, this effect can become moderately important 
and enhance the neutrino fluence in the very high energies 
[l[ . But this effect is hardly important in the case of the 
RS emission, as seen in Fig. 2. 

In Fig. 3, we show resulting spectra in the RS model 
with the ISM. For all the models, we hardly expect neu- 
trino signals by IceCube. As an example, in the model 
ISM-tc, we have small photomeson production efficiency 
f PJ ~ 5x 10~ 3 (at e p w 10 8 5 GeV). In the model ISM-tn, 



ivy 

smaller E Q j and smaller f pl 



5x10"* (at e p w 10 9 GeV) 



where e v is the neutrino energy in the comoving frame 
and g?£ is the luminosity distance to the source for a 



lead to much smaller neutrino fluences. In the original 
analytical prediction [27J , the neutrino energy fluence in- 
creases as neutrino energy is high. But, the more realistic 
spectra show the suppressed behavior around 10 9 GeV. 
This is just because protons have the finite maximum en- 
ergy so that we see the cooling effect of pions and muons 
in the highest energies. While the fluence of highest en- 
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FIG. 2: Neutrino fluences from one energetic GRB event at 
z — 1.6. Neutrinos are produced by the decay of pion and 
muon whose origins are single- pion production (thick lines), 
and double- and multi-pion production (thin lines). The 
model ISM-e is used. We show fluences of muon-neutrinos 
(i>n from pion decay (dot-dashed lines), muon-neutrinos 

{vp + P M ) from muon decay (dashed-lines) and electron- 
neutrinos (v e + i} e ) from muon decay (dotted lines). The total 
neutrino fluence is also shown by the thick solid line. The 
normalization of the proton flux is given by e aC c = 1/4 and 
neutrino oscillation is not taken into account. 



FIG. 4: The same as Fig. 3, but for the RS models with the 
wind-like CBM. 
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FIG. 3: Muon-neutrino (f M + £> M ) fluences (thick lines) and 
electron-neutrino (i>e + v e ) fluences (thin lines) from one GRB 
event at z = 0.1. The RS models with the ISM are used. Note 
that the normalization of the proton flux is given by e acc = 1/4 
and neutrino oscillation is not taken into account. 



ergy neutrinos is suppressed above E v ~ 10 9 GeV com- 
pared to the analytical result, we can see E 2 ^> v oc E],- 2 
for E v < 10 9 GeV, which is good agreement with the 
analytical result. In the model ISM-eb, we can see that 
neutrino spectra are extended to higher energies. This 
is just because the maximum energy of protons in this 
model, where large e r B = 0.2 is taken, is higher than in 
the models ISM-tc and ISM-tn (by a factor of ~ 5 larger 
than in the model ISM-tc). 
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FIG. 5: The same as Fig. 3, but model parameters are chosen, 
motivated by the modified RS model explaining the shallow 
decay emission. 



In the model ISM-tc2, we consider overlapping of the 
prompt emission. Once this is taken into account, neu- 
trino fluence is greatly enhanced by two or three orders 
of magnitude. Therefore, the effect of this overlapping is 
important. Nevertheless, the detection by IceCube will 
be very difficult unless a burst with E c j <~ ifj 54-55 ergs 
occurs at <~ 10 Mpc. 

As shown in Fig. 4, we can expect that the neutrino 
emission occurs efficiently in the model WIND-tc, where 
we obtain f P1 > 1 (for e p > 10 8 GeV). It is because the 
crossing radius in the model WIND-tc is smaller than 
that in the models with the ISM. Below E v - 10 8 GeV, 
we can see El4> u oc E^ 2 (note that the break correspond- 
ing to e™ w 8.3 eV is E™ - 10 9 8 GeV, which we cannot 
see because the fluence is saturated at energy such that 
/p7 ~ 1). While the spectrum of the model WIND-tc is 
in the fast cooling regime, that of the model WIND-tn is 
in the slow cooling regime. We cannot expect high pho- 
tomeson production efficiency in the latter model. We 
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FIG. 6: The same as Fig. 3, but for the models based on the 
late prompt emission model. The normalization of the proton 
flux is given by £ acc = 10. 



have / P7 ~ 3 x 1(T 4 (at e p « 10 85 GeV) due to the 
smaller E e - } and larger crossing radius. 

If we consider overlapping of the prompt emission for 
the models with the wind-like CBM, neutrino fluence can 
be higher by one order of magnitude due to additional 
target photons from the prompt emission. 

In Fig. 5, neutrino spectra are calculated under the 
RS model with the higher tB and smaller f e . We choose 
such parameters based on the modified RS model that 
explains the shallow decay emission [5l|, H3] ■ We calculate 
for the RS produced by the head of the ejecta, although 
the RS emission would continue by the slower tail of the 
ejecta. Although additional protons supplied by the tail 
of the ejecta may produce neutrinos, our estimation on 
neutrino fluences would not be changed so much up to a 
factor. Target photon spectra in both models ISM-s and 
WIND-s, are expected in the fast cooling regime at the 
crossing time of the head of the ejecta. In this RS model, 
significant energy of the RS emission is radiated as x- 
ray photons that can interact with protons with energy 
E p ~ 10 9 GeV, which produce neutrinos with £„~5x 
10 GeV. For the model ISM-s, we have the photomeson 
efficiency f pl ~ 0.01 (at e p « 10 85 GeV), which is higher 
than that in the model ISM-tc. It is because photons 
that can interact with protons increase in this RS model 
due to its higher injection energy ~ 0.3 keV. On the 
other hand, in the model WIND-s, we obtain f pl ~ 0.1 
(at e p w 10 9 GeV), which is smaller than that in the 
model WIND-tc. This is because most of the radiated 
energy is emitted as photons with ~ 7 keV which is 
higher than that in the original RS model. After all, it 
reduces the number of target photons that can interact 
with protons. Corresponding to e^b ~ ? keV, we expect 
the break at E™ ~ 10 7 GeV, which can be seen in Fig. 5. 
We can also find that overlapping of the prompt emission 
enhances neutrino fluence greatly in Fig. 5. 

In Fig. 6, we show results obtained under the late 



FIG. 7: The muon-neutrino (v^+v^) background from early 
afterglows of GRBs, based on the original reverse-forward 
shock model which typically predicts infrared/optical flashes. 
The normalization of the proton flux is given by e acc = 1/4 
and neutrino oscillation is taken into account. 

prompt emission model. For all the three models, we ob- 
tain fjyy > 1. Expected muon events from one GRB event 
is Nfj_ ~ (0.3 — 0.6) events for the model LP2. Therefore, 
if bright GRBs occur, we have possibilities to detect neu- 
trino signals from early afterglows in the late prompt 
emission model. Because most of the sufficiently high 
energy protons are depleted in these three models, we 
see the similar level of the energy fluence. In our mod- 
els, the magnetic field strength is stronger than that in 
the RS model. Hence, neutrino spectra are suppressed 
in the highest energies due to cooling processes of pions 
and muons. (Note that the model LPO has the same pa- 
rameter set used in Murase & Nagataki [53| , but we have 
the higher energy fluence by a factor because we showed 
the energy fluence per internal collision in the previous 
Letter.) 

B. GRB Neutrino Background 

Because it is difficult to see neutrino signals from one 
GRB event, we have to observe many GRBs. Because we 
can find a good fraction of GRBs by space satellites, we 
expect neutrino signals that are coincident with bursts 
seen by electromagnetic observations. Such a correla- 
tion is important to detect neutrino signals from GRBs. 
Otherwise, neutrino signals may be buried below the at- 
mospheric neutrino background and cosmogenic neutrino 
background in the high energies. 

We calculate the neutrino background from GRBs for 
our specific parameter sets. 

In Fig. 7, we show the neutrino background for the 
parameter sets based on the original RS model which 
was developed in the pre-Swift era. Even in the model 
WIND-tc, expected muon numbers by IceCube are ~ 
(0.05 — 0.1) events/yr. It is much difficult to see signals 
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FIG. 8: The same as Fig. 7. But predictions are based on 
the late prompt emission model and reverse-forward shock 
model with the higher es and smaller / e , motivated by Swift 
observations. The normalization of the proton flux is given 
by £acc = 10 for the late prompt emission model, while by 
e a cc = 1/4 for the reverse-forward shock model. 
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FIG. 9: The same as Figs. 7. But the RS model where 
the possible overlapping of the prompt emission is taken into 
account. The normalization of the proton flux is given by 

Cacc = 1/4. 



for the other models in Fig. 7. Even worse, if a RS is not 
common, neutrino signals are also even less expected. 

In Fig. 8, we show the neutrino background based on 
the late prompt emission model and RS model where 
the higher es and smaller f e are assumed. Both are 
developed for explanation of the shallow decay emis- 
sion. In the late prompt emission model, we can expect 
that the optical depth for photomeson production is high 
enough. In such cases, although UHECR production is 
difficult, GRBs are efficient neutrino emitters. We ob- 
tain ~ (0.5 — 1) events/year for the model LP1 and 
N/j, ~ (1 — 3) events/year for the model LP2, respectively. 
Although such signals can marginally be detected by Ice- 
Cube, the significant detection requires neutrino detec- 



FIG. 10: The optimistic predictions of the neutrino {y e + 
v e + + + v T + i? T ) background from early afterglows 
of GRBs in the Swift era. The upper three curves are for 
the model LP2, middle three curves for the model WIND- 
s2, and lower three curves for the model ISM-s2. The de- 
pendence of the neutrino backgrounds on GRB rate mod- 
els is also shown. Solid curves are for the GRB3 model, 
dashed curves for the GRB2 model, and dotted curves for 
the GRB1 model (see Appendix B). In the late prompt emis- 
sion model LP2, the normalization of the proton flux is given 
by Cacc = 50. In the RS models ISM-s2 and WIND-s2, 
the possible overlapping effect due to the prompt emission 
is taken into account, and the proton flux is normalized to 
E%dN p /dE p (z = 0) = 0.5 x 10 44 ergsMpc" 3 yr" 1 which corre- 
sponds to the observed UHECR flux. 



tors larger than IceCube. For the RS model, we have 
Nf, ~ (0.002 - 0.005) events/year for the model WIND-s. 

In Fig. 9, we show the neutrino background in the 
RS models where the possible overlapping of the prompt 
emission is taken into account. The effect of overlapping 
enhances neutrino fluxes greatly below ~ 100 PeV. 

In Fig. 10, we adopt the different normalization of 
the proton flux from Figs. 7-9. For the RS models, it 
comes from the hypothesis that the observed UHECRs 
come from the RS. Note that large baryon load is as- 
sumed for all the curves in Fig. 10. In addition, we have 
also assumed that all the GRBs are in thick ejecta regime 
for the RS models. Although such assumptions may be 
optimistic, they can lead to detectable neutrino signals 
from early afterglows by future observations. The depen- 
dence of the neutrino backgrounds on GRB rate models 
is also shown. We can see that the resulting neutrino 
backgrounds are not senstive to adopted GRB rate mod- 
els GRB1-GRB4 (where curves for the model GRB4 are 
not shown, but they are also similar to the other curves) . 

In Figs. 11 and 12, we show various predictions on the 
neutrino background from GRBs, motivated by observa- 
tions in the Swift era. Fig. 11 represents our fiducial pre- 
dictions. For cosmological (high luminosity (HL)) GRBs, 
we can expect ~ (6 — 11) event s/yr, which is consis- 
tent with previous works 0, IE] • For low luminosity (LL) 
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FIG. 11: The fiducial predictions of the neutrino (u e +9 e +u fl + 
Vp+Vt + Vt) background from GRBs in the Swift era. Prompt 
(HL): prompt neutrino bursts from cosmological (high lumi- 
nosity) GRBs; Bgrb = 1.24 x 10 51 ergs, E™° h = 10 51 ergs, 
r = 10 13 " 14 - 5 cm and T = 10 2 ' 5 , £s = 1 and £ acc = 10 0,0. 
Prompt (LL): prompt neutrino bursts from low luminosity 
GRBs; the local observed rate Pll(z = 0) = 500 Gpc _3 yr _1 , 
ET ~ 10 50 ergs, r = 10 15 cm, r = 10, £s = 1 and £ acc = 10 
[54| . Late Prompt: flaring neutrino flashes and/or neu- 
trino early afterglows under the late prompt emission model; 
E LP = O.ISgrb, L b = lO^ergss" 1 , r = 10 15 ' 3 cm, To = 15, 
£ B = 1 and Cacc = 10 (the model LP0) [Hj]. AG (WIND): neu- 
trino early afterglows under the reverse-forward shock model 
with the wind-like CBM; the model WIND-s is assumed (see 
the text). AG (ISM): neutrino early afterglows under the 
reverse-forward shock model with the ISM; the model ISM- 
s is assumed (see text). WB: Waxman-Bahcall bounds [77j . 
£ s and £ acc are defined as £ B = Ub/U~, and £ acc = U p /U 1 , 
respectively. For the fast cooling case and the proton acceler- 
ation efficiency C, p ~ 1, we have £s ~ (es/e e ) and £ acc ~ l/fie, 
where e e is the fraction of the internal energy density carried 
by electrons. In all the scenarios, the GRB3 model is used 
(see Appendix B). 



GRBs [5J, UM, we can expect N„ ~ (2 - 5) events/yr, 
although neutrino signals coincident with LL GRBs are 
usually not expected by Swift and GLAST. For flares 
[5oT ] and early afterglows under the late prompt emission 
model (the model LP0), we can expect ~ (0.4 — 1) 
events/yr. In the very high energy region of the neutrino 
background, neutrinos from LL GRBs, flares and early af- 
terglows will be dominant to those from HL GRBs. For 
AG (WIND) (the model WIND-s) and AG (ISM) (the 
model ISM-s), we expect ~ (0.002 - 0.005) events/yr 
and ~ (0.0002 - 0.0004) events/yr, respectively. 

Fig. 12 expresses the optimistic predictions. For HL 
GRBs, we show the neutrino background for £ acc = 50. 
The optimistic model for £ acc = 100 is being con- 
strained by pre sent and future-coming neutrino obser- 
vations [80 . l8l| . Note that the neutrino flux is al- 
lowed to exceed the WB bound in principle if the col- 
lision radii are small enough, but the neutrino flux for 
Prompt (HL) would typically not exceed the WB bound 



FIG. 12: The same as Fig. 11, but the optimistic pre- 
dictions are shown. £ acc = 50 is used for Prompt (HL) 
and Late Prompt. For Prompt (LL), AG (WIND) and AG 
(ISM), the proton flux is normalized to EpdN p /dE p (z = 0) = 
0.5 x 10 44 ergs Mpc _3 yr _1 which corresponds to the observed 
UHECRflux. 



[J [73, LZ8[- For Prompt (LL), we normalize the proton 
flux to E%dN p /dE p (z = 0) = 0.5 x 10 44 ergsMpc"V -1 
just for comparison. Although the shown parameter set 
allows protons to be accelerated up to the ultra high en- 
ergy fi p ~ 4 x 10 20 eV, the explanation of the observed 
UHECRs might be difficult if UHECRs are all protons 
[54| . (But note that heavier nuclei could be accelerated 
up to the highest energies.) If accelerated protons can- 
not achieve the ultra high energies (below ~ 10 19 eV), 
the neutrino flux does not have to satisfy the WB bound 
in principle and could exceed the flux shown in Fig. 12 
if the larger baryon load or higher rate is possible. For 
Late Prompt, UHECR production is typically impossi- 
ble, and we just adopt £ a cc = 50. For the RS models, we 
normalize the proton flux to the observed UHECR flux. 
We have - (0.2 - 0.5) events/yr for AG (WIND) and 
Nfj, - (0.02 - 0.05) events/yr for AG (ISM). 



V. SUMMARY AND DISCUSSION 

In this paper, we have studied the neutrino emission 
from early afterglows of GRBs under the assumption that 
baryon acceleration occurs at the shocks such as the RS 
and late internal shocks. With acceptable parameters, we 
have shown that neutrino signals from early afterglows 
could marginally be detected by IceCube under the late 
prompt emission model, while they are not expected un- 
der the RS model. Hence, the neutrino detection from 
early afterglows will be likely to suggest that neutrinos 
come from the late prompt emission (including flares). 
Future neutrino telescopes larger than IceCube are ben- 
eficial, although the simultaneous electromagnetic multi- 
wavelength observations are indispensable. 

Our conclusions are summarized below. 
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(1) We have revisited the neutrino emission from the 
RS. We have also taken into account the cross section 
of photomeson production quantitatively without using 
the A-resonance approximation. The effect of pion- 
multiplicity and proton-inelasticity is not important in 
th RS model and the A-resonance approximation is a 
good approximation. The neutrino flux is suppressed in 
the highest energies because of the finite proton's max- 
imum energy and cooling of secondary particles. In the 
original RS model, it is very difficult to detect neutrino 
signals for our fiducial parameter sets, even when all the 
GRBs accompany the RS. In addition, recent observa- 
tions imply the lack of infrared/optical flashes. If it is 
attributed to the intrinsically weak RS, we can expect 
neutrino production in the RS region only for a fraction 
of GRBs with the sufficiently strong RS. 

(2) One of the suppression mechanisms for in- 
frared/optical flashes is the overlapping of the prompt 
emission. In this scenario, the RS occurs, but cooling of 
electrons suppresses the emission in the infrared/optical 
bands. We have also taken into account this effect, which 
can dramatically enhance neutrino fluxes. In the case of 
the thick ejecta colliding into the ISM, enhancement of 
the flux by about two or three orders of magnitude is ex- 
pected, while by one order of magnitude in the case of the 
thick ejecta colliding into the wind-like CBM. Expected 
muon events can increase by a factor. Although this com- 
ponent may be hidden by prompt neutrino bursts, such 
contribution will be important in the sense that we do 
not have to assume the proton acceleration in the region 
where the prompt emission occurs. 

(3) Recently, it is suggested that the shallow decay 
emission in early afterglows may be explained by the 
modified RS model. As such an example, we have esti- 
mated the neutrino flux under the RS model with f e < 1. 
If this picture is true, we expect that the RS is com- 
mon for all the GRBs. Therefore, we can cumulate 
the neutrino flux from each burst and obtain the neu- 
trino background. We expect muon events by IceCube, 
Nfj, ~ (0.002 - 0.005) events/yr in the model WIND-s 
with the moderate baryon load. The hypothesis that the 
observed UHECRs come from the RS leads to the higher 
neutrino flux, but it requires the large baryon load. 

(4) One of the recently suggested models is the late 
prompt emission model. The physical conditions are sim- 
ilar to those in flares. With acceptable parameters, neu- 
trino production can be efficient because the emission oc- 
curs at smaller radii in this model. We can have / P7 > 1, 
which could allow us to estimate the proton flux from the 
neutrino flux. In this model, we could expect detectable 
muon events by IceCube, although the larger telescopes 
would be desirable for sufficiently significant detections. 
If we see such signals, neutrinos will provide us with use- 
ful information on baryon acceleration in GRBs and one 
of the clues to the model of early afterglows. 

(5) From Figs. 11 and 12, we have seen that the con- 
tribution from Prompt (HL) to the neutrino background 
may be the most important below ~ 10 PeV, while the 



neutrino background from early afterglows can be more 
important in the very high energies > 10 PeV. Neutrino 
signals from early afterglows are expected to be coinci- 
dent with the early afterglow emission and give us addi- 
tional chances to detect high energy neutrinos. However, 
note that the photomeson production efficiency / P7 is 
sensitive to collision radii r 0. For Prompt (HL) with 
To > 100, we have typically min[/ P7 , 1] ~ (0.1 — 1), and 
the neutrino background can be smaller than that shown 
in Figs. 11 or 12 where we have / P7 ~ 0.5 effectively (see, 
e.g., the curves for the parameter set B used in Murase 
& Nagataki [l[). On the other hand, we obtain typically 
/p7 ~ 1 f° r small T ~ 10, which is expected in case 
of the late prompt emission model with £"lp ~ 0.1-Egrb- 
Therefore, the neutrino background from early afterglows 
in the late prompt emission model can be comparable to 
that from the prompt emission. Roughly speaking, the 
three possibilities (Prompt (HL), Prompt (LL) and Late 
Prompt) can give comparable contributions to the diffuse 
neutrino background. 

It is important to see neutrinos associated with 
gamma-rays. For HL GRBs, flares and early afterglows, 
we can expect conincidence, while not for LL GRBs ex- 
cept for very nearby events. For the short-lived RS mod- 
els where the overlapping of the prompt emission occurs, 
it will be difficult to distinguish between RS neutrinos 
and prompt neutrinos. This is because the short-lived 
RS emission for the thick ejecta case lasts for the dura- 
tion of bursts ~ T. If the long-lived RS emission occurs, 
the long-term neutrino emission will be also expected. 
However, it may be contaminated by neutrino flashes 
from flares or the long-term neutrino emission from the 
FS. The latter could occur if protons are accelerated up 
to very high energies by other mechanisms such as the 
second-order Fermi acceleration (2^, [3(1 [M| . 

Our predictions are for specific parameter sets. Of 
course, parameters such as E b , E™° have dispersion. A 
possible wide range variability may also suggest disper- 
sion of E™° h . Some of the parameters might be related 
like the E b -E™° relation for the prompt emission. More 
comprehensive studies are needed in future, and more 
and more observations would enable us to take into ac- 
count the distribution of parameters. 

So far, the amount of protons has not been well con- 
strained from observations. In Fig. 11, we have shown 
the relatively moderate cases, while the optimistic cases 
are shown in Fig. 12. Due to the lack of knowledge on 
the collisionless shocks and particle acceleration, we can- 
not say whether such optimistic choices are possible or 
not. But note that too large values might be implausi- 
ble from recent observations [82]. Unless there is signif- 
icant missing energy (which might exist, for example, if 
accelerated protons carry away significant energy), high 
radiative efficiency of GRBs implies e e ~ (0.1 — 1), which 
leads to £ acc < 10. On the other hand, the UHECR hy- 
pothesis and recently estimated local GRB rate require 
the large baryon load for both of the prompt and RS sce- 
narios. We have assumed that nonthermal protons are 
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contained, which is expected from the internal or exter- 
nal shock model. However, this might not be true if the 
outflow is Poynting-dominated. The prompt or RS emis- 
sion may come from magnetic dissipation processes such 
as reconnection, where much baryons do not have to be 
contained. Future neutrino observations are fruitful, be- 
cause they could provide us with not only information on 
the physics of GRB outflows but also one of the clues to 
unknown baryon load and acceleration in GRBs. 

In this paper, we have focused on high energy neutri- 
nos from GRBs. Generally, the high gamma-ray emis- 
sion from proton-induced components such as muon syn- 
chrotron and secondary pair synchrotron should accom- 
pany neutrinos. Such high energy gamma-rays can ap- 
pear after complicated pair-photon cascades. Theoret- 
ically, efficient proton acceleration might induce dis- 
tinctive GeV-TeV components and such possibilities are 
studied by several authors [H, HH, [85|, [86|, [87| • If the over- 
lapping of the prompt emission is common, we would 
expect GeV-TeV flashes from the RS via the electron 
inverse- Compton scattering process [50| . In addition, the 
high energy emission via this process is predicted under 
the several early afterglow models motivated by Swift ob- 
servations H, H3| . Future-coming GLAST observations 
[9(il ] would also be useful for testing the models. 
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APPENDIX A: METHOD OF CALCULATION 

We briefly describe our method of calculation. The 
method is essentially the same as that used in Murase & 
Nagataki pHH. 

The most promising acceleration process of radiating 
particles is the Fermi acceleration mechanism. Electrons 
and protons are usually accelerated in GRBs by the Fermi 
acceleration mechanism. The first-order Fermi acceler- 
ation occurs by diffusive scattering of particles across 
strong shocks. For non-relativistic shock acceleration, 
the acceleration time scale is, 

3r 1 

— — (k u + r c Kd), (Al) 



where r c is the compression ratio, U\ is the velocity of the 
upstream fluid and n Ul d is the diffusion coefficients. In 
the ultra-relativistic shock limit, with the assumption on 
the sufficient amplification of the downstream magnetic 
field, we can obtain the acceleration time as (75| . 
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where r rc i is the relative Lorentz factor between fluids, 
i co h is the coherent length of the upstream magnetic field 
and tl is the Larmor radius. Let us assume the Bohm 
diffusion coefficient and sufficiently large coherent length 
of the upstream magnetic field. Then, we can write the 
acceleration time of protons as follows, 
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For the case of mildly relativistic (r re ] ~ a few) shocks, 
77 ~ 1 may be possible. The acceleration time scale for 
the second-order Fermi acceleration is given by, 
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where Pa = Va/c, va is Alfven velocity and SB is the 
strength of the turbulence in the magnetic field [73|, [9l[ . 
The limit 3a = 1 and SB = B which could be expected 
in the downstream of relativistic shocks also leads to the 
similar expression to Eq. (A3). 

Proton's maximal energy is also constrained by var- 
ious cooling processes. In this paper, we treat proton 
synchrotron cooling, inverse Compton (IC) cooling, adi- 
abatic cooling and photohadronic cooling. First, the syn- 
chrotron loss time scale for relativistic protons is, 

Second, the IC cooling time scale is given by [92| . 
272 \m 
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Third, the photohadronic cooling time scale is, 
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where s is the photon energy in the rest frame of pro- 
ton, 7 p is the proton's Lorentz factor, K,p IS the proton- 
inelasticity, and £th is the threshold photon energy for the 
photohadronic process in the rest frame of the incident 
proton. In sufficiently high energies, photomeson cooling 
process is much important, where the threshold energy is 
e t h ~ 145 MeV. We calculate Eq. (A7) by using Geant4 
whose total cross section is fairly good agreement with 
experimental data [55[ , and the calculated mean free path 
is alsogood agreement with that obtained by the other 
code [93, 94]. Fourth, we take into account the adiabatic 
cooling process, which has a time scale i a d independent 
of the proton energy. In fact, direct ejection of protons 
from the emission region may depend on a proton energy 
if diffusive losses are relevant. For simplicity, we neglect 
such diffusive losses and assume that protons are confined 
over the time scale set by adiabatic expansion. When the 
fluid is relativistic, we have, 



a ( 1 
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where V is the fluid velocity. From above time scales, 
the total proton loss time scale is expressed as t~ 1 = 
+ t~y n + tjQ 1 + . The proton's maximum energy 
can be determined by t acc < t p . 

Accelerated protons interact with target photons via 
photomeson production. Pion spectra can be obtained 
by executing Geant4, which are written as, 
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where dn p /de p and dn/de is proton and photon distribu- 
tion in the comoving frame, £„• is the pion-multiplicity. 
Owing to Geant4, we can include the effects of proton- 
inelasticity and pion-multiplicity. Although the treat- 
ments of photomeson production are greatly improved 
by using Gcant4, compared to the A-resonance approx- 
imation which is often used by other authors, Geant4 
has some problems in parametrization [l[ . Therefore, we 
adopt the more improved approximate treatment [531.154 
Below 3 GeV, we use the experimental data 5f| 
the previous paper [l|, we approximated ir + : tt° 
for single-pion production, and ir + Tr~ : ir + ir° = 7:4 for 
double-pion production pA [73[. Both approximations 
do not change our conclusions). Above 3 GeV, we use 
the Geant4 approximation. Such an approximation is 
sufficient for astrophysical applications to calculation of 
GRB neutrinos, and our results for neutrino spectra are 
quantitatively improved compared to most of the previ- 
ous works, where the A-resonance approximation is used. 

We can obtain neutrino spectra from the well-known 
decay kinematics. Neutrinos are produced by the decay 
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7m T M' respectively. Here, 7V = 2.6033 x 10~ 8 s and r M 
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where, e* — 



ilarly, we can get the muon spectrum from the pion 
spectrum. Muon decay is the three-body-decay process, 
which is slightly more complicated than the case of two- 
body-decay. Given the spectrum of muon, it can be cal- 
culated by the following equation [95j, 
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where, e* x = 7 M e y - (7? - 1) £„, e y2 



ill ~ 1) 1/2 ^, {ml - ml)c 2 /2m^ e* v - — 
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are defined in the muon rest frame and and for muon 
neutrinos, f {x) = 2a; 2 (3 - 2x), fi{x) = 2x 2 (l - 2x) 
and for electron neutrinos, fo{x) = 12a: 2 (1 — x), fi(x) = 
12x 2 {l — x), where x = 2e*/m AI c 2 and 9* is the angle 
between the muon spin and the direction of a neutrino. 

Because of cooling processes of ir^ and fi^ , we have to 
apply Eqs. (AlO) and (All) at each time step, when we 
solve the following equation, 
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where Q^,^ represents the source term of pions and 
muons due to photomeson production and decay of pi- 
ons, respectively. The synchrotron cooling time scale is 
given by replacing proton mass with pion or muon mass 
in Eq. (A6). The adiabatic cooling time scale is still 
comparable to dynamical time scale. We also treat the 
IC process including the Klein-Nishina effect for pions 
and muons. We neglect the 7r7 process (including the 
photomeson production and photopair production), be- 
cause other cooling processes are usually more important 
in our interested cases. Throughout the calculations, we 
also neglect neutrinos from neutron decay n — > p+e~ +9 e , 
whose time scale is usually much larger than the dynam- 
ical time scale idyn for the short-lived emission (but note 
that these components will also contribute to the diffuse 
neutrino background). 



APPENDIX B: GRB NEUTRINO BACKGROUND 

In order to get the differential number flux of back- 
ground neutrinos, first we compute the present number 
density of the background neutrinos per unit energy from 
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one GRB. The contribution of neutrinos emitted in the 
interval of the redshift z ~ z + dz is given as, 

dnf(E v ) = R GRB (z)(l + zf^-dz dNv ^ v) dE v {\ + z)" 3 , 



(Bl) 
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1 



where 

dt _ 

Tz - "mi + z) y /n A + n k {i + z) a + n m {i + z) 3 ' 

(B2) 

and E v = (1 + z)E v is the energy of neutrinos at red- 
shift z, which is now observed as E v and dN v (E v )/dE v 
is the number spectrum of neutrinos emitted by one GRB 
event. Hence, the GRB neutrino background can be cal- 
culated using the following equation, 



cdn° h c f z — ^ dN l/ ((l + z)E v ) 
- ' dzR GRB (z) 



dE v dVl AnH J 



dE' 



^n A + n k (i + z) 2 + a m {i + z) 3 



(B3) 



In this paper, we adopt z roax = H, &m = 0.3, ilk = 0, 
il A = 0.7 and H = 71kms _1 Mpc" 1 . Here, R GRB (z) 
is the beaming-corrected (overall) GRB rate. The lo- 
cal GRB rate has some uncertainties depending on the 
rate history, but our results on the neutrino background 
are not so changed because the main contribution to the 
background comes from GRBs that occur at z ~ (1 — 3), 
the number of which is observationally determined. We 
use the following GRB rate history in units of yr _1 Gpc -3 

QUi, 

46 e 3Az 

Rgrbi(z) = l8———F(z,n m ,n A ) (B4a) 
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R G rb 2 {z) = 18 „ . ■ F{z,n m ,n A ) (B4b) 
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94 „3.05z-0.4 

i?GRB 3 (^) = 23 2 My ; - F(z,Sl m ,Q A ) (B4c) 
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Rgrba(z) = 43F(z,Q m ,n A ) ^ ^0.75 [ z > 1} ( B4d ) 



where F(z,fi m ,Q A ) = yfljv + ^m(l + + ^) 3/2 - 

In this paper, we show results calculated under the GRB3 



model in most figures, because the dependence of the 
resulting neutrino backgrounds on GRB rate models is 
not so large [l[. In the previous work [l|, we evalu- 
ated the neutrino background by exploiting Eqs. (B4) 
with the beaming-corrected energy -Egrb = 1-24 x 10 51 
ergs. Alternatively, when we have isotropic energy, we 
should use the apparent GRB rate which is expressed as 
Pgrb = fbRcKB- In this paper, we adopt f b = 1/75 as a 
fiducial value for evaluation of the neutrino background 
from the RS. 

APPENDIX C: VACUUM NEUTRINO 
OSCILLATION 

Neutrino physics has been interesting since neutrino 
experiments recently discovered something new, rather 
than giving only more precise measurements of standard 
model parameters, or stronger bounds on unseen new 
physics. Solar and atmospheric neutrino data directly 
show that the flux of neutrinos with each flavor is not con- 
served, which suggests that neutrinos are massive then 
there should be a neutrino mixing matrix. 

In GRBs, neutrinos are produced via photomeson pro- 
duction, so that we can expect v e : : v T sa 1 : 2 : 
approximately. As a result of neutrino oscillation, we can 
obtain v e : : v r w 1 : 1 : 1 . So there may be a possibil- 
ity that tau neutrinos are detected through double bang 
events 97]. However, we should note that the magnetic 
field is strong in GRBs, so that contributions from muons 
that are more subject to cooling due to their longer life 
time, are suppressed. Hence, we expect v e : : v r w : 
1 : in the high energy region. As a result of neutrino 
oscillation, we can obtain v e : : v T « 1 : 1.8 : 1.8 
rather than v e : : v r sa 1 : 1 : 1 in the high energy 
region (9ST |. 

In this paper, we consider vacuum neutrino oscillation 
in the long baseline limit. Assuming that 023 is maximal 
(6*23 ~ 7r/4) and #13 is very small (#13 w 0), which are 
indicated by neutrino oscillation data, we can obtain, 

*u. « <-isin 2 20 12 (2< c -$^~O (CI) 

^ « *„ r « 5(^ + 

+ i sin 2 2M2< - - < ). (C2) 
In this paper, we adopt #12 = 0.59. 
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